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Abstract 

For a SU {N) Yang- Mills theory, we present variational calculations using 
gaussian wave functionals combined with an approximate projection on gauge 
invariant states. The projection amounts to correcting the energy of the 
gaussian states by substracting the spurious energy associated with gauge ro- 
tations. Based on this improved energy functional, we perform variational 
calculations of the interaction energy in the presence of external electric and 
magnetic fields. We verify that the ultraviolet behaviour of our approxima- 
tion scheme is consistent, as it should, with that expected from perturbation 
theory. In particular, we recover in this variational framework the standard 
one-loop beta function, with a transparent interpretation of the screening and 
anti-screening contributions. 
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I. INTRODUCTION 

The functional Schrodinger picture has proven to be a priviledged tool in exploring a 
rich variety of aspects of gauge theories which are beyond the scope of perturbation theory 
It is a useful starting point for developping non perturbative calculations based on 
the variational approach. In the case of scalar field theories, static as well as dynamical 
variational calculations have been performed by using trial wave functionals of the gaussian 
type In the case of gauge theories, some early investigations along these lines can 

be found in [pHT^. However the application of variational methods to gauge theories is 
generally plagued by the difficulty to implement in a calculable way the requirement of 
gauge invariance of physical states |T^-|T3[. 

Gaussian wave functionals allow for analytic calculations, but are not gauge invariant 
except in the Abelian case. In principle, one can construct gauge invariant states by averag- 
ing gaussian wave functionals over all gauge rotations. This results in an effective non-linear 



sigma model where the fields are the group elements of the gauge transformations How- 
ever, to make progress with this theory, further approximations are necessary both in the 
choice of the kernel of the gaussian, and in the evaluation of the functional integrals over the 



gauge group [p!^ -[18|. Such approximations, which go beyond the variational principle, are 
not always under control. In particular, in the perturbative regime, they fail to completely 
reproduce the one- loop beta function [jl5| (see also |[T6| , p!8| ). 



In this paper we shall propose a different strategy which is inspired by techniques devel- 



oped in 1962 by Thouless and Valatin [|T9[ to deal with the restoration of rotational invariance 
when deformed solutions are obtained in nuclear Hartree-Fock calculations. Rather than us- 
ing gauge invariant variational states, we shall limit ourselves to gaussian wave functionals, 
but we shall correct the associated energy functional by a non-local term, which approxi- 
mately corresponds to the energy gain when projecting on gauge invariant states. In the 
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Abelian case, this amounts to removing the contribution of the longitudinal part of the 
gaussian kernel to the energy. In the Yang-Mills theory, the corrective energy term is itself 
determined by the variational principle, and the ensuing variational calculation is a priori 
non-perturbative. 

Our main purpose here is to provide the framework for such a calculation and verify 
that, in the perturbative regime (i.e., when supplemented with an expansion in powers of 
g), it leads to results consistent with ordinary perturbation theory. As a first check, we 
shall show that our variational calculation reproduces the standard one-loop beta function 
of Yang-Mills theory, with a transparent interpretation of the various contributions in terms 
of screening and anti-screening phenomena. As a further test, we shall compute the vacuum 
energy in the presence of an external magnetic field B, and find that it exhibits a minimum 



at a non vanishing value of B, in agreement with the perturbative calculation in pO . 

The paper is organized as follows. In Sec. II we briefly review the functional Schrodinger 
picture and the variational principle in field theory. In Sec. Ill, we present the Thouless- 
Valatin formalism and consider, as a simple illustration, its application to quantum elec- 
trodynamics (QED). In Sec. IV we present a variational calculation of the one- loop beta 
function which is based on the construction of the interaction energy between external elec- 
trostatic charges. In Sec. V we consider the energy of the QCD vacuum in the presence of 
a constant magnetic field. This provides an alternative computation of the beta function, 
and also of the gluon condensate which is found to satisfy the trace anomaly relation. A 
summary of our results and a discussion of further possible extensions and applications is 
presented in Sec. VI. 

II. VARIATIONAL CALCULATIONS FOR GAUGE FIELDS 

We consider the functional Schrodinger description of the SU{N) Yang-Mills theory. 
In the temporal gauge = 0, the canonical coordinates are the vector potentials v4^(x) 
and the electric fields ii^*(x), which we shall often write as color matrices in the adjoint 
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representation: e.g., = A\T^ (the color indices a, b, . . . run from 1 to A^^ — 1). The 
generators of the color group are taken to be Hermitian and traceless; they satisfy 

The Hamiltonian density reads {g denotes the coupling constant) : 

7^(x) = i i^g'ElEli^) + ^5:i?:(x)| , (II. 1) 

with the color magnetic field: 

B„ = V X A« + ^ fabcA^ X Ae. (II.2) 

Note that our conventions are such that the QCD coupling constant is absorbed in the 
normalization of the vector potentials. With these conventions, the covariant derivative 
reads Di = di — iAi, and the electric fields i?* are canonically conjugate to the vector 
potentials A^: [E^(x), Ai{y)] = i6'^6ab6'^'\^ - y). 

In the Schrodinger representation, the states are represented by functionals of A^(x), 
\1/[A], and the electric field is acting on such states by functional differentiation: 

= .^vE.[A]. (II.3) 

The Hamiltonian H commutes with the generator Q of time- independent gauge transforma- 
tions, 

^?(x) = V-E(x) + t[A,E% (II.4) 

so it is possible to diagonalize H and Q simultaneously. 
The physical states are constrained by Gauss' law: 

gi^)^[A] = 0, (II.5) 

which is the requirement of gauge invariance. [More precisely, eq. ( [11. 51 ) shows that physical 
states must be invariant under "small" (i.e., topologically trivial [Q) gauge transformations. 



We shall not be concerned with the topological aspects of the gauge symmetry in what 
follows.] More generally, in the presence of matter fields represented by an external color 
source with density p", Gauss' law gets modified as follows: 

Gai^mA] = Pai^MA]. (II.6) 

The ground state of QCD is the eigenstate "^vac of H of minimal energy which satisfies 
Gauss' constraint ( [II. 5|) . It can be constructed, at least in principle, by using the Ritz 
variational principle, which states that: 

with the minimum achieved for \1/ = \l/^ac- Here, "^[A] is any wave functional from the 
physical Hilbert space (i.e., which satisfies eq. ( [II.5|) ), and E^ac is the energy of the ground 
state "^vac, assumed to be non-degenerate. In practice, however, one has to restrict oneself 
to gaussian wave functionals, the only ones which allow an analytical computation of (if). 
These have the form 

^o[A] = A^-iexp{-^/d3xd3|/[Aj,(x)-4(xy 

(11.8) 

where the background field ^^(x) and the kernel G^^ (with matrix elements (G^^)"* (x, y)) 
are the variational parameters. We expect that A = in the vacuum state, and this is the 
case that we shall consider mostly in this paper. Still, non-vanishing values of A will be 
also considered, in Sec. V below, in a study of the vacuum stability in the presence of a 
background color magnetic field (in the same spirit as, e.g., in Ref. ||20|| ) . 



The expectation value of the Hamiltonian density in the gaussian state \l/o is |T2| 
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(^ol^(x)l^o) = 2^B-B(x) + |Tr(x|G-i|x) 
+ |Tr(x|ii'G'|x) 

(II.9) 

+ ^ (Tr{S,r'^(x|G|x)}f 

+ ^Tr {^^T'^(x|G|x)^^T'^(x|G|x)}. 

In this equation B is the magnetic field associated to the center A and S** is the spin one 
matrix whose elements {j, k) are given by ieijk- The notation Tr in equation ( p.9|) implies 
a summation over the discrete (color and spatial) indices. For instance, Tr(x|G'^^|x) = 
J2i,a (^~^)i'r(x5 x). Finally, the operator K is the second derivative of the classical energy 
with respect to the center A";. It reads, in matrix notations, 

K = {-tS-T>f -S-B, (11.10) 

where 

V' = d' - iA' (11.11) 

denotes the covariant derivative defined by the background field A^ = A^T"". In particular, 

Kij{p) = p^5ij - PiPj for A = 0. (11.12) 

In the case of non-Abelian gauge theories, however, gaussian functionals like eq. ( [1I.8|) 
suffer from a major drawback: they do not satisfy the requirement of gauge invariance ( |1I.5| ). 
It is in principle possible to construct gauge invariant states by projection, i.e. by averaging 
a gaussian functional over all its gauge transformations. This is achieved by means of the 
formula 

vI/[A] = -l|l?[f/(x)]vI/^[A], (11.13) 
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where the functional integration is performed over the unitary N x N matrix field f/(x), 
with the adequate group invariant measure, and A/" is a normalization factor. The integrand 
in eq. ( [11.13| ) is the gauge-transform of the gaussian state ^^olA]: 

^t/[A] = ^o[UAU+ + iUV (11.14) 

The expectation value Ep of the energy in the projected state is given by the following 
formula 

_ JP|t/(x)|<1-o|g|$„> 
JVlU(x)]<%\^tu> • 



which should replace eq. ( [1I.9|) in practical calculations. Unfortunately, eq. (|1I.15|) cannot 



be evaluated in closed form because the functional integral over the group is not gaussian. 
Various approximations to eq. (p.l5|) have been considered in Ref . |[T^-[T8| . In what follows, 
we shall propose a different approximation method which is inspired from techniques used in 
nuclear physics to calculate the zero point rotational energy of deformed nuclei [p!9| , |2T| -p^ . 

III. APPROXIMATE PROJECTION 



The fact that our variational ground state, namely the gaussian \&o in eq. (|II.8|) , is not 



gauge invariant introduces a spurious degeneracy in the problem: \E'o is degenerate with all 
its gauge-transforms defined in eq. ( |1I.14| ). This leads to the existence of spurious excitations 



of zero energy which corresponds to rotations of \E'o in the gauge space; these are, of course, 
the Goldstone bosons associated to the spontaneous breaking of the gauge symmetry by 
\E'o. Accordingly, the expectation value (p.9|) of the Hamiltonian in the deformed state \I'o 
includes unphysical contributions expressing the kinetic energy of the gauge rotations. The 
Thouless-Valatin formalism |19[ provides us with a method to estimate, and thus subtract 
away, such unphysical contributions. 



7 



A. The case of rotations 



This formalism is best explained on the example of the collective rotations of a deformed 
nucleus. There, the equivalent of our present variational calculation with gaussian wave 
functionals is the so-called Hartree-Fock approximation where the nuclear wave function is 
represented by a Slater determinant formed with A single-particle wave functions (pk{x) (for 
A nucleons) The latter are determined by solving the Hartree-Fock equations, i.e., the 
variational equations obtained by minimizing the expectation value of the Hamiltonian in the 
subspace of Slater determinants. Although the Hamiltonian H is rotationally invariant, the 
Hartree-Fock equations may lead to non-invariant solutions describing nuclear deformations. 
If this is the case, then the spectrum develops a ground state rotational band, 

Ej = Eo + ^^^^, (II1.1) 

where J{J + 1) is the eigenvalue of the angular momentum operator and / is the moment 
of inertia. One may then conclude that the operator: 

H^H-^ (111.2) 

is the Hamiltonian of the intrinsic (that is, non-rotational) motion: to the ground state 
rotational band of H corresponds now a single eigenvalue Eq of H. Then, eq. ( [111.2| ) provides 
an approximate separation of the dynamics into intrinsic and rotational motion, which is 
reminiscent of the familiar separation of the center-of-mass motion by the formula: 

H = H , (III.3) 

where 

P = Pi + P2 + . . . + PA, (III.4) 

is the total momentum operator, and M = Am is the nuclear mass. Thus, a mean field 
description of the intrinsic motion can be given by performing variational (or Hartree-Fock) 
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calculations for the substracted Hamiltonian H. The only question is, what is the value of 
the moment-of-inertia parameter / ? 

One can answer this question by studying rotations of the deformed mass distribution. 
Assume that the nucleus has axial symmetry with respect to the z-axis, and consider an 
uniform rotation with angular velocity ui about the x-axis. If ^uj{t) is the exact (time- 
dependent) state describing such a rotation, then 

^^(t) = e-''"*'^-e-^^"**^(0), (III. 5) 
which satisfies the time-dependent Schrodinger equation, 

= H^^, (III.6) 
provided ^a)(0) is a solution to the following time- independent problem: 

{H - ^ J,)*a;(0) = ^.^..(0). (III.7) 

This is, of course, just the familiar transformation to the rotating frame of reference, which 
leads us to consider the variational problem for the following Hamiltonian: 

H^^H - ujJ^. (III.8) 

This is equivalent to a constrained variational problem with the subsidiary condition that 
< Jx> has a given value; uj plays then the role of the Lagrange multiplier. The constrained 
Hartree-Fock calculation determines the optimal independent nucleon wave function '^^ = 
^u>{^) in the rotating frame. (As a; — > 0, — > which is the deformed Hartree-Fock 
ground state.) Once the optimal wave function is found it is possible to obtain an estimate 
of the moment of inertia / of the nucleus by considering the limit 

I = li^SlAM^^_ (IIL9) 

This finally allows us to estimate the zero-point rotational energy AEtv in the deformed 
state '■ 
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Eq. ( [111.101) is the expected gain in energy when projecting the deformed Hartree-Fock 
ground state \l/o onto a rotationally invariant state. That is, the corrected average energy 
after projection, which is the energy of the intrinsic motion (cf. eq. ( [I11.2| )), reads 

E = Eo-AEtv, (III.ll) 

where Eq =< ^ol-f^l^o >• 

The Thouless-Valatin method is an approximation which is expected to be vahd for large 
deformations, or more precisely when the deformation produces a large expectation value of 
the square of the angular momentum in the Hartree-Fock ground state Indeed in this 
case it can be shown that the projection onto invariant states can be accurately performed 
because the overlap between two states differing in their orientation by an angle 6 is sharply 
peaked near 6=0 thus allowing for an an expansion in the vicinity of this point p2| . 



B. Application to QED 

When going to gauge theories in the variational method, the Hartree-Fock ground state 
\E'o is replaced by the gaussian trial wave functional (cf. eq. ( [II. 8|) ), and the angular momen- 
tum operator is replaced by the generator of the gauge transformations, ^"(x) (cf. eq. (p.4[ )). 

As a first illustration let us consider a variational calculation in QED, with the gaussian 
variational Ansatz 

vl/o[A] =Ar-'exp{- <A\^\A >}, (III.12) 

where J\f = {detCy^^ and the expression in the exponent is a condensed notation for the 
convolution 

<A\^\A>^y d'xd'yA\^)Gr^\^,y)A^{y). (III.13) 

This wave functional is gauge invariant provided its kernel G^^ is transverse: diG^j^ = 0. Let 
us assume, however, that this is not the case, and see what the Thouless-Valatin correction 
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would predict in this case. The operator playing the role of the angular momentum is the 
charge density operator 

g(x) = V-E(x), (III.14) 

and the generalization of the Thouless-Valatin formula ( [1II.10| ) for the energy correction 
reads 

AEtv = J d^xd^y < ^o\G{^)Q{y)\'^o > < x|^|y >, (III.15) 

where the "moment of inertia" X is now a matrix in coordinate space, with matrix ele- 
ments < x|X|y >= X(x, y). This is obtained via a constrained variational calculation with 
Hamiltonian H^^ = H — Hext, where 

H = ^J&'x {E2(x) + (V X A)2}, (III.16) 

and the external constraint 

Hext = J d3x^(x)V-E(x). (III.17) 

In the present context, the Lagrange multiplier co'(x) plays the role of the temporal com- 
ponent A'^(x) of the gauge vector potential. The solution to this constrained variational 
problem is of the formQ 

*,[A] = exp{-i < F|A >}^o[A], (III.18) 

where the vector field F(x) is a new variational parameter, which expresses the expectation 
value of the electric field in the state ( [111.18| ), F* =< '^^{E^l'^i^ >, and is determined by 
minimizing 



"^Eq. (III. 18) can be simply understood by recalling that, in the presence of a constraint of the 
form Hext = ctp, the ground state wavefunction of an harmonic oscillator is modified by a factor 
exp {ipQx), but its width remains unchanged. 
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E^ = <<iI^\H- H,^t\-^^ > (III. 19) 



2J ' '2 
We have denoted here {V is the total volume of the space) 

Eq = < ^ol^l^o > (III.20) 
V f d'p ri 



2 



/ (2^ [I'^ii^iP) + -m)G'u(p)} 



Note that the magnetic piece of the energy ( |1II.20|) (the second term between parentheses) 
involves only the transverse components of Gij, while the electric piece involves also its 
longitudinal component. 

The functional -^^[F] in eq. ( [1II.19|) attains its minimum for F = — Vcu, in which case 

< ^^|V-E|^^ > = V-F = -Acj. (111.21) 

According to ( |111.9|) , the moment of inertia is obtained as (with < V-E >^ = < 
^^|V-E|$^ >) : 

< x|J|y > = - 



5uj{y) 

= <x|-A|y>, (III.22) 
whose inverse is simply the Coulomb propagator: 

<x|y|y >= -. (III.23) 

I 47r|x — y| 

We thus obtain the following expression for the Thouless-Valatin correction ( [111.15| ) in QED: 

AEtv = \j d^xd^^y 

< ^ol V-E(x) V-E(y)|^o > — ^ ; , (111.24) 

47r|x — y| 



which is recognized as the electrostatic energy in the state \t'o- For the gaussian state (pi.l2| ), 
this gives 
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which simply subtracts the longitudinal piece of the electric energy in eq. ( |III.20|) . 



The corrected energy E = Eq — AEtv reads therefore (cf. eqs. ([111.2CI| ) and ([111.25| )) 



^ = ^ / (W [l^^'^P^ + p'GAp)} , (111.26) 



d^p f 1 

and involves only the transverse piece Gt = ^{^ij — PiPj)Gij of the kernel Gij (we have 
written here p = |p| and pi = Pi/p)- Then, the variational equation SE/SGt = gives Gt 



in the expected form 



Grip) = ^ . (III.27) 



Together, eqs. ( [ill.2(j| ) and ( [111.27] ) yield an energy density 



E^VJ^^M, (III.28) 

which is indeed the correct result for the QED ground state Thus, in the case of QED, the 
approximate projection method of Thouless and Valatin correctly subtracts the contribution 
of the unphysical, gauge, degrees of freedom from the average energy. Actually, since QED 
without fermions is a free theory, the variational solution above coincides with the exact 
solution 1^: the exact ground state is a gaussian wave functional like eq. ([111.12| ) with a 
transverse kernel determined by eq. ( 1II.27|) : 



S.aJAl = exp |- I ^ I A'{p)(S„ - p,p,)A'i-P)] 

^exp{-/d'.d-v|™). (III.29) 

Moreover, wave functionals of the type shown in eq. (|111.18| ) — i.e., gaussian states with a 
transverse kernel and a non-trivial phase factor — correspond to physical charged states, i.e., 
states of the quantum Maxwell theory in the presence of static, classical, external sources. 
Indeed, any such a state (which we denote here as \E'c) satisfies: 

V-E(x)^,[A] = p(x)^,[A], (111.30) 

with the charge density p(x) = V-F(x). The corresponding energy includes the Coulomb 
energy, as expected: 

13 



= E„ + i/d'xd'y|^^W^. (111.31) 

2 J 47r|x — y| 

Non-Abelian charged states will be considered in Sec. IV.B below. 

For other applications and a more complete study of the Thouless-Valatin method in the 
context of quantum field theory, see p4|. 



C. Approximate projection in QCD 

Let us now consider the case of non-Abelian gauge theories. The corresponding "moment 
of inertia" is now a color matrix defined as the polarization tensor (cf. eq. (|II.4| )) 

X^^(x,y) = ^^|J|^Uo, (III.32) 
in the presence of an external constraint 

Huj = H — Hext 

Hext = J d^x cu»(x)6?'^(x). (III.33) 

The analog of eq. ( [111.1(J| ), i.e. the gain in energy when projecting a wave functional \l/o[A] 
onto the subspace of gauge invariant states, reads: 

AEtv = I d'xd'y{<ilo\G''i^)g\y)\^o){a, x|^|6, y) (III.34) 

The energy functional to be used in variational calculations is therefore 

E = Eo- AEtv (III.35) 

It is also possible to perform a projection on a subspace with a given color charge dis- 
tribution < ^"(x) >c7^ (The subscript c refers to expectation values over charged states). 
In such a case, however, the gauge generator appearing in the previous correction formula 
has to be replaced by its deviation Q away from the desired value: 
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^-^(x) = g-{^)- < 6;«(x) >, . (III.36) 

This modification guarantees that there is no correction for a state which is an exact eigen- 
state of the charge operator. The functional to be minimized in the subspace of gaussian 
functionals is thus 

E = E,- Eext - AEtv, (111.37) 

where E = {H)c, -Eext = (-f^cxt)c is the energy of the external constraint generating the 
charged stateQ, and 

AEtv = I d^xd'y < g^i^W'iy) >c (a,x|^|6,y). (III.38) 

This procedure is again reminiscent of the elimination of the center-of-mass motion in the 
mean field description of a composite system of A particles For a system characterized 
by a set of single-particle wave functions (pi, (p2, ^Pa the optimal state in the center-of- 
mass frame is obtained by minimizing the functional 

E={m-^, (111.39) 

The Thouless-Valatin prescription for the total mass M in eq. ( [111.39| ) is to use the relation 
(P) = Mv, where (P) is the expectation value of the total momentum ( |111.4| ) in the presence 
of the external constraint Hext = v • P. This prescription gives the desired result M = mA 
where m is the mass of the individual constituents. In a moving frame with velocity v the 
single particle wave functions become 

VJi(x)->e*>^(^V^(x), (III.40) 

with = ■ X- Individual momentum operators in the moving frame are obtained by 
the gauge transformation 



^That is, an exact charged state is defined as an eigenstate oi H — Hg^f 
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V + 2(Vx), (111.41) 

i.e. Pi Pi — (pi). The functional providing the adequate state at the minimum is 

E={H)-Y-P- ^ (III.42) 

in agreement with eq. (|III.37|) . This procedure to implement Gauss's law will be important 
in Sec. IV when applied to the calculation of the interaction energy of color charges. 

Still in the case of charged states, the chromostatic energy Echromo is given, in our ap- 
proximation scheme, by the classical chromostatic energy corrected by the Thouless-Valatin 
term : 

Echromo = J d'xd'y{g''i^))c{g\y))c{a,^\^\b,y) + AEtv (111.43) 

so that 

Echromo = J d'xd'y {g%^)g\y))c{a,^\^\b,y). (III.44) 

(This is the analog of using Erot =< > /2/ as an approximation for the rotational energy 
of a deformed nucleus.) In the case of QED this identification is obvious on equations like 

( Pt:21 or ( pmp . 



To conclude, the central result of this section is the non-local energy functional (|III.35|) 
(or ([III. 37] ) in the case of charged states) which approximately corrects for the lack of gauge 
symmetry when working with gaussian states. This energy functional is the starting point 
of the variational method we propose for gauge theories. Note that the corrective energy 
term AEtv in eq. ( [III.34|) is a priori of a non-perturbative nature. Our aim in what follows 
is to check the ultraviolet behaviour of this approximation scheme. We shall thus consider 
the variational calculations in the perturbative regime g <^1. 



IV. ONE-LOOP BETA FUNCTION FROM VARIATIONAL CALCULATIONS 

In this section, we shall use eq. ( [1I1.43| ) to estimate the electrostatic energy Echromo of 
a non-Abelian charged state, up to order in perturbation theory. This will allow us to 
recover the standard expression for the QCD beta function in the one-loop approximation. 
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A. Moment of inertia for color rotations 



The first step is the calculation of the moment of inertia for color rotations, Tafe(x, y), 
which enters eq. (p.ll.43| ). Unlike QED, where this quantity has been computed exactly (cf. 
eq. ( |111.23|) ), in QCD we shall give only a perturbative estimate of /, valid to the order of 
interest (i.e., up to order g"^). To this aim, it is sufficient to perform variational calculations 
in the vicinity of the perturbative vacuum. 

To zeroth order, the vacuum of the Yang-Mills theory is the same as for the Maxwell 
theory, namely (cf. eq. ( [111.291 )): 



^o[A] = AT-i exp <^ - ( A 



A 



(IV.l) 



where 



(G-^)^^(k) = 5., U, - ^) G-^\ 



(IV.2) 



and G^^ = 2k. In the calculation of / below, we shall never need to go beyond this leading 
order approximation for Gk- 

Note that, even with such a transverse kernel, the functional ( [IV. 1|) is still not invariant 
under non-Abelian gauge tranformations; that is, this is a deformed state, according to the 
terminology in Sec. III. In order to compute its moment of inertia under color rotations, 
one has to study the response of this state to an external constraint of the form ([111.331) . 
The trial wave functional in the presence of this constraint reads : 



Ag' 



A 



(IV.3) 



which involves two additional variational parameters: the vector field F^{x.) (which fixes the 
expectation value of the electric field: (E^) = F^), and the matrix S^^^(x, y), which is taken 
to be transverse in its spatial indices. The emergence of S is a hallmark of the non-Abelian 
behaviour (recall that S = in the corresponding Abelian problem; cf. eq. ( p.ll.l8[ )). 
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We shall shortly see that, as — > 0, F and S are linear in c<j, while the first cor- 
rection to is only quadratic^, and therefore does not matter for the calculation of X 
(cf. eq. ( [111.32| )). Thus, for the present purposes, we can take G^^ as in eq. ( [1V.2| ). The 
variational parameters in ( [1V.3| ) are obtained by minimizing the following functional 

= <H J d^xuj'' <g'' (IV.4) 
with respect to variations in and E^J,. A straightforward calculation yields: 

J d^x{F, ■ F,(x) + VTr(x|SGS|x)}, 



(IV.5) 

where we have kept only the terms involving the variational parameters. Similarly, 

< ^"(x) = V • F"(x) - i/Tr(x|T''[S, G]\x.). (IV.6) 



After inserting these expressions in eq. ( |1V.4|) , and taking variations with respect to F and 



S, we derive the following expressions for the variational parameters (in momentum space) 



i^^(q) = -^^a(q), (IV. 7) 



and, for the transversePI components of S, 



{k' b\J:\k c) = ^^5{k' -k-q) (IV.8) 



-'- a/ \ r / Gk — Gk' \ 



^Such quadratic corrections occur since, in contrast to QED, the non-Abelian constraint ( pi.33| ) 
generates a coupling between G and F in the variational equations. 
'^By which we mean transversality with respect to both k and k', as requested by the expressions 



(FVTSI) and (FV6|) 



where V is the total volume and is defined in (p!V.2| ). 

By using eqs. ( [1V.6| ), ( [IV .7] ) and ( [IV .81 ), we can finally express the average color charge 
< Q >uj in terms of cj". This is conveniently decomposed into an "Abelian" and a "non- 
Abelian" piece, as corresponding to the two pieces in the r.h.s. of eq. ( [1V.6| ): < Q"" >uj= 
Pa + Pna^ with: 

1 



and, respectively. 



p^(q) = • F'^(q) = - q2a;'^(q), (IV.9) 

9^ 



PiV^(q) = 9 fabc J (5,, - -^)(% - -^) 

X (Gk - Gk') < k'c|S|k6 >, (IV.IO) 



where k' = k + q. From eqs. ( P^V.8[) -( [[V.10|) we note that the "non- Abelian" charge density 
in eq. ( [IV. 10 ) is a correction of order g"^ relative to the "Abelian" one in eq. ( IV. 9 ). Thus, as 



anticipated after eq. ( [IV. 2|) , it is consistent to evaluate this correction with the free kernel 
Gk = l/2k. 

The integral in eq. ( |1V.10 ) is logarithmically ultraviolet divergent, so it must be evaluated 



with an upper cutoff. It turns out that this divergence is a part of the charge renormalization 
in QCD (see Sec. IV.B below). To reconstruct the associated beta function, we need, as 
usual, only the coefficient of the divergent logarithm. The latter is insensitive to the details 
of the UV regularization, so we shall consider, for simplicity, a sharp momentum cutoff A. 

Also, in order to isolate the leading logarithm, we can perform kinematical approxima- 
tions relying on the inequality k ^ q (since the external momentum q is fixed, while the 
leading contribution to the integral in eq. (P^V.10[) comes from relatively large momenta). 
Physically, we are indeed interested in smooth charge distributions. This allows us to re- 
place k + q by k and thus {6ij — ^)(5jj p^) by 2. By also using eq. ( p.V.8[ ) for S, we 



then obtain 
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where Cn = Tr (T^T^) = N for SU{N), and 

X(a) - f [^(kO - ^(k)]^ ,2) 

= J (27r)3e(k).(k')(.(k)+.(kO) ' ^^^-^^^ 

with e(k) = |k|. Here again we can replace k' ~ k everywhere except in the numerator 
which must be expanded to second order in q: 

[£(k + q) -£(k)]2 -g^pos^^. (IV.13) 

where 6 is the angle between the space vectors k and q. The angular average yields (cos^ 6) = 
desmOcos^e = 2/3, so, finally: 

(IV.14) 

which is logarithmically divergent in the ultraviolet, as expected, but also in the infrared: the 
infrared divergence is an artifact of the previous manipulations (in a more careful calculation, 
this would be screened by g), and to the order of interest we can just regulate it with an 
ad-hoc infrared cutoff /i. This yields = ap% with 



a^^ln^. (IV.15) 

487r2 /i2 V ^ 



The resulting value of the moment of inertia lab = 6p'^/6uj reads finally: 

Xa6(q) = %g'(l + a). (IV.16) 
9 

This should be compared to the corresponding Abelian result^ I{q) = q^. We see that 
the quantum fiuctuations in QCD produce an increase of the moment of inertia, which 
corresponds to the screening of color charges by quantum fiuctuations. The size of the 
screening effect that we have obtained agrees indeed with the results of other approaches 
2|. 



^The overall factor in eq. ( [IV.16| ) is simply a consequence of our different normalizations for 
the field strengths in QCD and QED; compare, in this respect, eqs. ( [III.12| ) and ( [V.l| ). 
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B. Interaction energy in the presence of a background electric field 



Interesting properties of the vacuum include its response to an external chromo-electric 
field, which can be generated by an external constraint of the form 



(IV. 17) 



For this constraint, we shall compute the induced electric mean field and charge density, and 
the associated electrostatic energy. By comparing the latter with the bare Coulomb inter- 
action, we shall then identify the chromo-electric susceptibility, or charge renormalization. 
As we shall see, the variational formalism provides a transparent picture of the underlying 
phenomena of screening and anti-screening. 

The optimal state \E'c in the presence of the constraint ( [iV.lTp is of the form 



1 „-i<F|A> 



exp 



A 



A 



(IV.18) 



where the parameter F"- (the electric mean field) will be related to the external field El^^ 
by the variational equations (see eq. ( [IV.23|) below). Note that, in contrast to eq. (|IV.3| ), 
there is no S term in eq. ( [IV. 18|) above; this is so because the external perturbation here 
is different (compare eqs. ( [IV.17|) and (|III.33|) ): it contains a term linear in E", but no 
non-Abelian term like This situation is analogous to the case of an anharmonic 

oscillator with an external constraint Hext = ctp. The only chang factor e*^"^ and, in 

higher orders in g, a modification of the real part of the width. 

According to the discussion in Sec. Ill (see especially eqs. ( |III.37| ) and ( p.II.38| )), the 
energy functional to be minimized in this case is E = Ec — (-f^ext)c — AEtv- The terms 
involving F^ in this functional read: 

1 



Ec — E, 



ext 



AE' 



TV — 9 



d^x <! ^ F" • F" 



"R" ■ F" 



x(x|G'^||y)(a,x|^|6,y) 



(IV. 19) 
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The last term in the r.h.s. corresponds to the Thouless-Valatin correction, eq. ( [III.38|) ). 
Note that, because of the substraction of the average charge in Q"" = Q"'— < >c, it is 
only the non-Abelian part of the Gauss operator (i.e., the term i[A\ E^] in eq. ( [11. 4| )) which 
contributes to eq. ( [1V.19| ). To evaluate this contribution, we first rewrite it in momentum 
space: 

d^A; d^g 



(27r)3 (27r)3 

To the order of interest, we can replace the moment of inertia in the equation above by its 
leading order expression: X(k) ^ k'^/g'^- Then, by performing similar manipulations as in 
the calculation of the X term in eqs. ( [1V.10D -( [1V.14D , we finally obtain: 

^^^^ " ^ y (2^0^ (2^0^ ^^^^^ ^-^^ 2k¥ 

jd\F^{^)F^{^\ (IV.21) 

where 



Gtt^ j k 127r2 y/i^ 
After inserting ( |1V.21| ) in ( [1V.19| ) and minimizing with respect to F°(x), we obtain 

(1-5)F^(x) = E:,,,(x). (IV.23) 

At this point it is convenient to introduce the charge distribution associated to the 
external field -Eexf 

P:,,(x) = V-E:,,(x), (IV.24) 

to be referred to as the external charge in what follows: this would be the charge in the 
system in the absence of polarization effects. The actual charge is rather 

= Y^^ pU^ + s)^ (IV.25) 
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where the second hne follows from eq. (|IV.24| ). Note that this relation implies an antiscreen- 



ing of the external charge, since p*^ is bigger than p^^^. The difference 

p — Pext = Pext S may be interpreted as an induced charge (see also Sec. IV. C in Ref. and 
Appendix A for an alternative computation of this quantity). 

We are finally in position to compute the chromostatic interaction Echromo in the optimal 
state \E'c- This is given by eq. (|III.43| ) which, together with the above expressions ( [1V.25|) 



for (^"(x))e, and ( WM for AEtv, implies: 

1 + 35 f d^g ^„ ^ 1 



Echromo ^ 9^ -—— / TT^-T^ Pext (q) Pext (-q) 

1 + a J [Zirr 



1 + a J (27r)3 ^^-t^^^^^-t^ 2q2 ' 

(IV.26) 

up to corrections of higher order in g. There is here a factor (1 + 5)^ ~ 1 + 25 arising from the 
induced charge (cf. eq. (|1V.25|) ), another one arising from the Thouless-Valatin correction 
( [IV. 211) , and a factor (1 + a) due to the moment of inertia (cf. eq. (|IV.16|) ). The interaction 



energy (|1V.26|) is still Coulomb like. 



EcHromo = qI j ^ Pext ( q) Pe^xt ( " q) ^ , (IV.27) 



but with a modified coupling constant given by 



or, to first order in g"^ 



aUf^) = (IV.28) 



' ' ^^^-m^. (IV.29) 



gUp) 9' 48vr2 p2 • 
This is the correct one-loop value for the renormalized coupling constant [^. Note that, in 
the present calculation, this involves three types of contributions: indeed, the factor 11 in 
the last equation has arised as 11 = 8 + 4— 1, where the 8 corresponds to anti-screening by 
the induced charge (cf. eq. ( [1V.25[ )), the 4 is another anti-screening contribution due to the 



Thouless-Valatin correction ( [1V.21| ), and the (—1) is a screening contribution arising via the 



correction of order g"^ to the moment of inertia. 

23 



V. VACUUM ENERGY IN A MAGNETIC BACKGROUND FIELD 



In the previous section, we have studied the electric sector of the vacuum of the Yang- 
Mills theory, by using a combinaison of variational and perturbative techniques. In what 
follows, we shall perform a similar analysis of the magnetic sector. To this aim, we consider 
the Yang-Mills theory in the presence of a (constant) magnetic background field B^, and 
compute the background field energy by using the variational principle. The final result is 
not new (it coincides with the one-loop result by Savvidy |^), but it rather serves as a test 
for our variational method in the magnetic sector, and in the perturbative regime. 

The relevant trial wave variational is the gaussian functional \l/o in eq. ( [II. 8| ) with the 
"center" field ^^(x) chosen so as to reproduce the desired magnetic field 5* (a convenient 
choice will be given later). This state is not gauge- invariant, so its energy Eq = (^o|-f^|^o)) 



eq. ( |II.9|) , must be corrected with the Thouless-Valatin energy AEtv^ to be computed in 
the next subsection. Then, by applying the variational principle to the corrected energy 
E = Eq — AEtv, we shall determine the kernel of the Gaussian (in Sec. V.B). Finally, 
in Sec. V.C, we shall compute the energy of the magnetic field and the associated gluon 
condensate, and verify that these quantities are related by the trace anomaly relation, as it 
should. In this calculation, the standard one-loop beta function will emerge once again. 



A. The Thouless-Valatin energy in the background field 



We start by computing the moment of inertia Xab(x, y) in the presence of the background 
field A^(x). As already explained, this requires constructing the variational ground state 



\E'a; for the constrained Hamiltonian H^^ in eq. (|III.33| ). This state is of the form (compare 
to eq. (irVl ) 



X exp 









4 







A- A 



(V.l) 
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where the parameters and Tjf- are related to uo"- by the variational principle (recall the 
discussion in Sec. IV. A). It turns out that the matrix S will not play any role in what 
follows: indeed, below we shall need the moment of inertia only to leading order in while 
S counts starting with order g"^ (cf. Sec. IV. A). We then write, as in eq. 



= <H>^-Jd^x c^"(x) < e;"(x) (V.2) 

with < H >^ given by eq. ( [IV. 5|) , and 

< 6;«(x) >^ = (P.FT (x) + 0(/), (V.3) 

where = 9* — iA^ is the covariant derivative defined by the background field (cf. 
eq. ( [11. Ill )), and the neglected terms, of 0{g^), would involve S (cf. eq. ( [lV.(j| )). The 
variation with respect to yields then 

^:(x) = --^{V^con^), (V.4) 

which differs from the corresponding expression in eq. (P^ V. 7]) only by the replacement of the 
ordinary derivative di by the covariant one Pj. Together with eq. ( [V.3[ ), this provides the 
moment of inertia to the order of interest: 

X"^(x,y) ^ ^yy!^ = A {Virs^'H^ - y) + Oil) 



6u'{y) g 



= ^(a, x|tf|6, y) + 0(l). (V.5) 
9 

We have introduced here the kinetic momentum Uj = iVj = i6j + Aj, and 11^ = Iljllj. 
Within the same accuracy, one has also: 

(vI/o|^"(x)^''(y)|vI/o) ^ ^,Vt:^^V^^ iG'')^fi^,y) 

= ^(a,x|n,Gr.in,|6,y). (V.6) 

We are now in position to compute the Thouless-Valatin energy AEtv, cf. eq. ([III.34|) : by 
combining eqs. ([V.5[) and ( [V.6[) , one obtains: 



AEtv - I J (a, x| (n.^H,) GT^'\a, x) , (V.7) 
up to corrections of order g"^. 
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B. The variational equation for G 



The improved energy functional E = Eq — AEtv reads therefore (cf. eqs. ( [11 .91 ) and 



E = Id'^l ^5,^^(x)Sf (x) 



+ ^ (a, x|(5,j - ni|^nj)G'i/|a, x) 

+ ^Tr[A-G(x,x)]+ (V.8) 



Note that the last two terms in ( [II. 9| ) do not contribute to this order. As obvious from this 
equation, the Thouless-Valatin correction makes the kinetical part of the energy covariantly 
transverse. Since the operator Kij is transverse as well (cf. eq. ([II.10[ )), 



Ki, = u% - UiU, + 2[n,,n,] , (V.9) 

it follows that the projected energy ( jV.Sj) involves only the transverse components of the 
kernel G. Thus, without loss of generality, we can restrict ourselves to a (covariantly) 
transverse kernel in what follows: 

rfG-.i = = G^'W. (V.IO) 

To formalize this, it is convenient to introduce transverse and longitudinal projectors as 
follows: 

4^n,ijn„ Q^i-P. (v.ii) 



They satisfy : 



11- 

P.2. = n — iii.n''' n-' = p 



Ql = {l-P)l = il-P),^=Q,^. (V.12) 

Then, a transverse kernel is one satisfying G = QGQ (and similarly for G^^). For such a 
kernel, the variational principle (i.e., the minimization of E, eq. ( [V.8[ ), with respect to G) 
produces the following gap equation: 
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which determines G to the order of interest. In particular, as A ^ 0, G reduces to the free, 
or Abehan, expression in eqs. (P^V.2|) and ( [III.27| ). Thus, the only non-trivial effects which 
are taken here into account are those associated with the background field. 

Note that eq. ([V.13| ) can only be valid at sufficiently high energy, or small coupling 



constant: indeed, the operator K admits negative modes p0|j27|j28|] . We thus assume that 



an infrared cutoff has been set - this does not affect the ultraviolet behaviour of the theory, 
which is our main interest here. 

C. The energy of the background field 

The previous equations provide the optimal gaussian kernel for a given background field 
A and thus the effective potential V{A) which is the expectation value of the energy in this 
state. The next step in our variational approach is to find the minimum of the effective 
potential. Constructing V for an arbitrary background is however a difficult task. For this 
reason we now consider a restricted variational space defined by the following background 
field : 

= ; Ay = xBT^ ; = 0. (V.14) 

This corresponds to a constant magnetic field in the z-direction and in the third color. 

With this choice of the background field, we are now able to compute the energy ( |V.8| ) 
in the optimal variational state, which is the gaussian state (|11.8|) with a transverse (in the 



sense of eq. ( |V.ICI| )) kernel G ^ satisfying eq. ( |V.13[ ). The latter equation shows that, at the 



minimum, the following identity holds: 

Tr{K(x|G'|x)} = ^Tr (xlG-^lx). (V.15) 

That is, magnetic and electric fluctuations have equal energies in our variational ground 
states, which is merely the virial theorem in the present context. Thus, 
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(V.16) 



This involves the matrix element (xjC ^|x), which we shall compute in Appendix A by using 
the Schwinger proper-time representation (cf. eq. ( |V.13| )) : 



1 roc rlf 

(x|G-i|x) = |^(x|(l-e-*^')|x) (V.17) 

This integral develops ultraviolet divergences as t 0, which we shall regularize by shifting 
the lower bounds of the integral from to As in the electric case, we are mainly 

interested in the ultraviolet renormalization of the energy ( |V.16| ); to this aim, it is sufficient 
to extract the terms which diverge when A — oo in eq. ( |V.17| ). This is described in detail 
in Appendix A, from which we quote here the final result: 

IY(x|G-^|x) = (. . .)A^ - ^ 5^ In ^ + 0{g% 

(V.18) 

(the coefficient in front of A^ is an uninteresting field-independent number that will be 
omitted in what follows). Remarkably, there is no divergent term in A^ (which, for di- 
mensional reasons, would be necessarily of the form K^B): this is so because of rotational 
and gauge symmetries which require the magnetic field to enter only in the scalar product 
S ■ BcjT" = BSzT^i whose trace is however zero (see Appendix A for more details). 

Note also the numerical factor in front of the logarithmic divergence in eq. (|V.18|) : this 
is the factor leading to the correct one-loop beta function after renormalization (see below). 
The projection on (covariantly) transverse gaussian states has been crucial in getting this 
factor right: without this, we would have obtained a factor | instead of the correct factor 



n 

3 



(compare in this respect eqs. ( |B.19| ) and ( |B.20| ) in Appendix A). 



Finally note that the field strength B appears as an infrared cutoff in eq. ( |V.18| ). This is 



expected from equations ( [B.6| ) and ( |B.7| ) where the proper time variable always appears in 
the combination tB. However a complete derivation of ( |V.18| ) requires a detailed treatment 
of unstable modes. 
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To conclude, 



E/V r:^ - CM^—\n—, (V.19) 



showing that the background field energy has no field dependent UV divergences other than 
the logarithmic one which can be absorbed into the renormalization of the coupling constant. 
We then write, as usual (/x is the substraction scale), 

11 11 

which provides the correct one-loop beta function, as anticipated. The renormalized field 
energy density reads then: 

which coincides with the result obtained by Sawidy in perturbation theory |^ . An advan- 
tage of the present approach, however, is that it can be improved by using a larger variational 
space, which is expected to cure the difficulties associated with negative modes P7|P^ . 
As discussed in ||2^ , the energy density ( [V.21|) exhibits a minimum for a non-zero value 



B = Bmin of the background field, with 



The value of the energy density at this minimum is : 

(^)mm = ^Q^^^^^rnin (V.23) 

which is indeed negative. Our variational vacuum state is therefore characterized by a 
magnetic field condensate (see, however, Refs. [0,^ for potential problems with such a 
state). 

From the previous results, it is now straightforward to evaluate the gluon condensate in 
our variational vacuum: 



1 

2i?Ln, (V.24) 
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where the second hne follows from the aforementioned "virial theorem" ([V.15|) . Eqs. ([V.23|) 
and ( |V.24| ) can be combined into: 

11 Cn 

which is consistent, as it should, with the trace anomaly relation: 

K) = ^ {F'-F,.). (V.26) 

Indeed, with (7i) = ^{O'^), and the one-loop beta function (which here is a consequence of 
eq. (^) 

m = -^,C^9' , (V.27) 

eq. ( |V.26D becomes identical to eq. ( |V.25| ). 

An attractive feature of the formula ( |V.25D is that it involves two quantities which are 
independently accessible experimentally (at least indirectly). Indeed, the left hand side of 
this equation is the energy density of the vacuum, which can be identified with the fourth 



power of the bag constant, B'^ = — (240Me\^)^ p9|, whereas the right hand side depends 
on the gluon condensate which is known from Ref. to be O.bGeV^. These values are 
compatible with eq. ( |V.25| ) within a 20 percent accuracy. 



VI. CONCLUSION 

In this paper we have proposed an improved energy functional for variational calculations 
in gauge field theories. This functional contains a non local term which approximately 
corresponds to the energy gain when projecting on gauge invariant states. This allows 
one to use gaussian states as trial functionals and thus perform analytic calculations for 
physical observables such as the chromoelectric and chromomagnetic susceptibilities, energy 
expectation values and the gluon condensate. 

The main purpose of this work was to check the ultraviolet behavior of our approximation 
scheme. By performing variational calculations near the perturbative vacuum we have shown 
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that divergences can be eliminated by a renormalization of the couphng constant. This has 
allowed us to recover the familiar one loop beta function in a way which makes transparent 
the various screening and antiscreening contributions. In particular the screening term arises 
naturally in our formalism, which was not the case in earlier variational approaches |]T5| . We 



have also tested our variational method in the magnetic sector, checking that it reproduces 
the one-loop result by Savvidy |^ for the background field energy. This calculation provides 



us with another derivation of the one-loop beta function. 

Thus, our formalism appears to correctly reproduce the expected behavior of non-Abelian 
theories in the ultraviolet sector. This strongly encourages us to study its predictions in the 
non perturbative regime. Indeed, as a variational approach, it is not at all restricted to 
the vicinity of the perturbative ground state, nor to small values of the coupling constant. 
We would like to also emphasize that the Thouless-Valatin correction is the first step in 
an approximation scheme which can be constructed systematically. Indeed, it is the first 



order term [g^ in boson expansion methods which have been constructed by Schwinger pT 
Dyson, Holstein-Primakoff Blaizot-Marshalek 



An attractive feature of the variational picture is that it allows one to treat situations 
where instabilities occur. This is the case when the trial state is not the lowest one and 
where saddle points are reached. An example of such a situation is the occurence of negative 



modes generally found in the presence of a constant magnetic background field ||27| , ^ . In 
this case it would be interesting to work out what is the optimal gaussian kernel in our 
approach. This state should be reached by allowing the variational space to include the 
subspace spanned by the negative modes. 

Another attractive question is the investigation of the infrared behaviour. It may provide 
some information on the confinement mechanism and the generation of mass scales . Sum 



rules and the gluon condensate at finite temperature also appear to be a promising field of 
investigation. 
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APPENDIX A: PHYSICAL INTERPRETATION OF SECTION IV.B. IN TERMS 

OF THE INDUCED CHARGE 

Let us show an other interpretation (and computation) of the relation between and 
the external charge. What follows is directly inspired from Gottfried and Weisskopf in 



section IV. C of |2^. Let us assume that the system is in presence of a given distribution 
of external charges plxti^) = ^ ' Eg^,((x) and compute the corresponding induced charge 
created by the quantum fluctuations of the gauge field ; it is given by the mean value < 
" fabcA'l{x.) {ElY (x.) > where E^, is the longitudinal part of the chromoelectric field operator 
E [^. This operator is a non-dynamical variable which is fixed by Gauss' law: 



V ■ E2(x) = pe.*(x) - UcAl{x){Em^) (A.l) 

which can be solved perturbatively, setting E2(x) = E2^°^(x) + E^^^^(x) + E^^^^(x) + .... In 
fact, only the first two terms will be needed for a development of the total charge in first 
order of g"^ . They verify the following set of equations : 

V-Er(x)=p:,,(x) 

V ■ E«'^(x) = -UcAlix)iElfH^). (A.2) 

The first equation shows that E^^(x) = Ee^t(x). Then the second equation is solved by: 

d^y Xi - yi 



(Eim^) = -J ^^^Jc,eAj{y){EL,ny) (A.3) 
The total charge reads therefore : 



= p:,,(x)+ < -/,,,4^(x)(i?Lj'=(x) > 

+ < -fa,cAl{x){ElfH^) > 

= P:,,(x)+p^„,(x). (A.4) 
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The first term in p^^^ vanishes since linear in Ai, while the second term yields 

pl,{^) ^ < ~UcAl{^)iElfH^) >, (A.5) 

or, after a Fourier transform, 

P-nd(q) = fa.JceiE^.tYM) J ^e'^^^ ^ G^/(x). (A.6) 

By also using 
we finally deduce: 

P-n.(q) = fa,cME^.t)-{c^) J ^ Glf{k). (A.8) 

Since G is transverse this expression reduces to 

Furthermore for a smooth charge distribution we can approximate k + q by k in the above 
integral. Noting that the result vanishes unless i=j we have 

PLd(q)=7P:.,(q), (A.IO) 

with 

1 rOD A h 



The total charge is then given by 



pL(q) = pL(q) + p"„d(q) 

= pUq){l + 6), (A.12) 



which is precisely the expression obtained in section IV. B (cf. eq. ( |IV.25|) ) 



33 



APPENDIX B: PROPER-TIME CALCULATION OF THE ENERGY DENSITY 



Let us present here in some detail the calculation of the quantity TrG"^ which enters 
the energy of the magnetic field in Sec. V.C According to the Schwinger proper-time 
representation ( |V.17| ), one needs the matrix element (x|e~*''^|x). Since, moreover, we are 
mainly interested in the ultraviolet behaviour of the energy, this expression is needed only at 
small values of t, which allows us to perform expansions in powers of t whenever necessary. 

As explained in Sec. V.C, we shall use the background field in eq. ([V. 14|) for which: 



[Vi, VjY^ = -pbs^sij-^ ^ T^S.B, (B.l) 

and therefore ^^,0^^'' = -[Vi,Vj]''^ = -{T'^)abiS,)ijB. It is convenient to define the 
following operator (cf. eq. (|V.9| )): 

kij = Kij + UiU, = U^6ij + 2[Ui, Uj] , (B.2) 

in terms of which Kij can be rewritten as follows : 

K = KQ. (B.3) 

By also using [K, Q] = 0, and Q + P = 1, we deduce 

We thus have to compute the matrix element {x\e~^^ Q\x.) , with : 

k = 11^ - 2T^S,B 
(x|e-*^|x) = (x|e-*n'e2*^'^^^|x), (B.5) 

where the second line follows since [n^,T^S'z] = 0. 

The computation of (x|e~*^^|x) is well-known in the literature with the result : 



Ant smh[tT-^B) 



We thus obtain: 
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(x|e-*^» = (^)3/2 tT'B 2tT35.B] _ 7) 

By also using Q = 1 — P and Hjll^ = KUi, one then rewrite (x|e~*^(5|x) as : 

(x|e-*^|x) - (x|e-*%in,i|x) (B.8) 

The last term of this equation can be obtained from (x|e~*''^njnj|x) by integration over t. 

To calculate (xle"*^!!^!!^ |x), we follow Schwinger's method : We work in Heisen- 
berg's representation with t = is and deduce 

(x|e-*"'nin,|x) = (xie-'^n^nin^ix) 

= (x(s)|n,(o)n,(o)|x(o)) (B.9) 

where : 

The operator U (s) = e"*"*^^ can be interpreted as the evolution operator of a particle gov- 
erned by the Hamiltonian 11^. We have 

= 2U{s) (B.ll) 

and 

^ = ^[U^U]{s) (B.12) 

Using [n2,n](s) = 2[nfc,ni]n'= = -2iJ^ikU\ where 

= i[Vi, Vjf^ = iT^S.B (B.13) 

one obtains n(s) = (e^'*-^)ifcn^(0), and thus : 

Xi{s) - x,(0) = ( ^ ),,n^(o) 

^R-'W{Q) (B.14) 
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The matrix element (x|e ^'^^Hjlljlx) can be now computed as (cf. eq. ( p.9|) ): 

(x(.)|n,(o)n,(o)|x(o)) 

= {x{s)\RimRjniXm{s) " (0) ) (x^ (s) - X„ (0) ) |x(0) ) 
= RimRjn{As)\Xm{s)Xn{s) - X„(s)x„(0) 

-Xn{s)Xm{0) +Xm{0)Xn{0) + [x„ (s) , (0)] |x(0) ) 
= i?,„i?,„i?-,l[nfc(0), X„(0)](X(S)|X(0)) 

= -2i?,,(x(s)|x(0)) (B.15) 
Returning to the variable s = —it, wee obtain 

(x|e-'"^n,n,|x) = ( ~f J .,(x|e-*"V), (B.16) 



and thus, finally, 



(x|e-*^ninj|x) 



^ ^ ^3/2 " ^tT^S.B 



2 sinh(tr3S^S) ' Ant ' smh{tT^B) 

-i^f^^l+T^BS^t 
2t^A-nt' \ 



+ - )t' + ■■■} (B.17) 

where in writing the second line we have developed up to the second order in t. 

The last equation can be now integrated over t to obtain the second term in (|B.8| ): 

(x|e~*^P|x) = {^y^'ll + T'BS^t (B.18) 

Developing (xje^^'^lx) up to the second order in t : 

(x|e-*^|x) = i-^fHl + 2T^BS,t (B.19) 
Airt [ 

+(2(r3)2525f - + . . . |, 
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we deduce the following expression for U{t) = (x|e *^Q|x) : 

f/(t) = (^p|^+T3B^.t (B.20) 



Thus, finally, 

IY(x|G-|x) = -^j ^(rr(x|g|x) -Trt/(t)) 



•)A'-^(^i?^)ln^ + 0(/), (B.21) 



3 ' B 



where (x|(5|x) has been obtained by setting t = in eq. ( |B.20| ). This is the result announced 
in eq. { ^^ ). 
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